home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
The X-Philes (2nd Revision)
/
The X-Philes Number 1 (1995).iso
/
xphiles
/
hp48_1
/
trap_dn
< prev
next >
Wrap
Internet Message Format
|
1995-03-31
|
3KB
From: David Nurser <itb625nurser@qut.edu.au>
Subject: v04i032: trap_dn - Trapezoidal Rule v1.0, Part01/01
Newsgroups: comp.sources.hp48
Followup-To: comp.sys.hp48
Approved: spell@seq.uncwil.edu
Checksum: 2883128221 (verify with brik -cv)
Submitted-by: David Nurser <itb625nurser@qut.edu.au>
Posting-number: Volume 4, Issue 32
Archive-name: trap_dn/part01
BEGIN_DOC trap.doc
Trapezoidal rule Program by david nurser (ITB625NURSER@redgum.qut.edu.au)
Works exactly the same as the built in integration (via stack)
(almost)
Accepts:
1: Lower Limit
2: Upper Limit
3: Equation (variable "X")
4: number of iterations
Produces:
Output...matrix of working (for exams)
h, the sum, and the answer
END_DOC
BEGIN_RPL trap.rpl
%%HP: T(3)A(D)F(.);
\<< 'X' PURGE \-> a b
EQ n
\<< b a - n / 0 \-> h
SUM
\<<
[[ 0 ]]
{ 'n+1' 3 } RDM 0 n
FOR r r 3 * 1
+ r PUT r 3 * 2 + a
r h * + DUP 'X' STO
PUT r 3 * 3 + EQ
\->NUM DUP
IF r 0 == r
n == OR
THEN 2 /
END SUM +
'SUM' STO PUT 1
STEP h "h"
\->TAG SWAP SUM "SUM"
\->TAG SWAP SUM h *
"ANS" \->TAG
\>>
\>> 'X' PURGE
\>>
END_RPL
BEGIN_ASC trap.asc
%%HP: T(3)A(D)F(.);
"D9D20E16324563284E20108597632EFE021C432D6E201016D6E201026D6E2020
5415D6E2010E6E1632D6E201026D6E20101690DA1D6E2010E650FA14B2A21C43
2D6E201086D6E20303555D4E16328E9209200033920200001000010000000000
000000000047A208BA20D6E2010E69C2A276BA1B21303F2A2B2130FD0D14B2A2
D6E2010E60A132D6E201027D6E2010273F2A2EEDA19C2A276BA1D6E201027704
D1D6E2010273F2A2EEDA1ED2A276BA1D6E201016D6E201027D6E201086EEDA17
6BA178BF14563284E20108597632DCC02704D1D6E2010273F2A2EEDA13F2A276
BA1D6E202054154E5A178BF13CE22D6E2010274B2A2279E1D6E201027D6E2010
E6279E1908E1AFE22D9D20ED2A250FA1B21305DF22D6E20303555D476BA14563
2D6E20303555D497632DCC02704D19C2A208332D6E201086C2A207000086EB52
2DBBF1D6E20303555D4C2A20B00003555D4EB522DBBF1D6E20303555D4D6E201
086EEDA1C2A20B000014E435EB522EF532EF5324563284E20108597632EFE029
3632B213047DF"
END_ASC
BYTES: #FD74h 397
BEGIN_UU trap.uue
begin 644 trap
M2%!(4#0X+46=+>!A(U0V@N0"`5AY-N+O(,$TTN8"`6%M+A`@UN8"`D51;2X0
MX.9A(VTN$"#6Y@(!80FMT>8"`6X%KT$K*L$TTN8"`6AM+C`P5=7D82/H*9`"
M`#,I(````0`0`````````````'0J@*L";2X0X)8L*F>KL1(#\Z*R$@/?T$$K
M*FTN$.`&&B-M+A`@U^8"`7+SHN+>&LFB<K8:;2X0('=`'6TN$"`W+RKNK>$M
M*F>KT>8"`6%M+A`@U^8"`6CNK7&V&H?[064C2"X0@)5G(\T,<D`=;2X0(#<O
M*NZM,2\J9ZO1Y@("15'DI7&X'\,NTN8"`7*THB*7'FTN$"#7Y@(!;G+ID8`>
M^B[2V0+>HE+P&BLQ4/TB;2XP,%75=+8:5#;2Y@(#4U5->3;2S"`'U)$L*H`S
MTN8"`6@L*G```&B^)=*['VTN,#!5U<2B`@L`,%75Y%LBO?O1Y@(#4U5-;2X0
E@.;>&BPJL```04Y3OB7B7R/^-4)E(T@N$("59R/^#I)C(RLQ`"LQ
`
end
END_UU